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■ ABSTRACT 

6 

Multi-object spectroscopic observations of blue supergiants in NGC 3621, a spiral 
galaxy at a distance of 6.7 Mpc, carried out with the ESO VLT and FORS are pre- 



sented. We demonstrate the feasibility of quantitative stellar spectroscopy at distances 
approaching a ten-fold increase over previous investigations by determining chemical 
composition, stellar parameters, reddening, extinction and wind properties of one of our 
targets, a supergiant of spectral type Al la located in the outskirts of NGC 3621. The 
metallicity (determined from iron group elements) is reduced by a factor of two relative 
to the sun in qualitative agreement with results from previous abundance studies based 
on H II region oxygen emission lines. Reddening and extinction are E(B — V) = 0.12 
and Ay = 0.37, respectively, mostly caused by the galactic foreground. Comparing stel- 
lar wind momentum and absolute V magnitude with galactic and M31 counterparts we 
confirm the potential of the wind momentum-luminosity relationship as an alternative 
tool to estimate extragalactic distances. 

Subject headings: galaxies: individual (NGC 3621) — galaxies: stellar content — stars: 
early-type — stars: winds, outflows 



1 Based on observations obtained at the ESO Very Large Telescope. 
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1. 



Introduction 



Detailed optical work with 4m-class ground-based telescopes and the Hubble Space Telescope 
on individual blue supergiants has been carried out in recent years in the Galaxy, the Magel- 
lanic Clouds and a handful of nearby galaxies, most notably M31 and M33 (McCarthy et al. 1995; 
Kudritzki 1998; Monteverde et al. 1997; Monteverde et al. 2000; Venn 1999; Venn et al. 2000a; 
Venn et al. 2000b; Dufton et al. 2000; Smartt et al. 2000). These stars are extremely luminous 
(up to My = —10) and show signatures of strong mass outflows in their ultraviolet, optical 
and infrared spectra, which are successfully interpreted in the framework of the radiation driven 
wind theory (Kudritzki & Puis 2000; Kudritzki 2000). One of the main predictions of the the- 
ory, nicely confirmed by spectral diagnostics of supergiants in the Milky Way (Puis et al. 1996; 
Kudritzki et al. 1999), Magellanic Clouds (Kudritzki 1998; Puis et al. 1996), NGC 6822 (Muschiel- 
ok et al. 1999), M31 (Smartt et al. 2000; McCarthy et al. 1997), and M33 (McCarthy et al. 1995), 
is a tight relationship between the modified wind momentum (the product of mass-loss rate, ter- 
minal velocity of the winds and square root of stellar radius) and the stellar luminosity: 



with a ~ 2/3 (Kudritzki 1998; Kudritzki & Puis 2000; Kudritzki 2000; Kudritzki et al. 1989). 
This Wind momentum - Luminosity Relationship (WLR) has a great potential as a new, inde- 
pendent extragalactic distance indicator. Kudritzki et al. (1999) estimate that using the WLR 
technique with ten to twenty A-supergiants per galaxy and medium resolution spectroscopy with 
8m telescopes from the ground distance moduli with an accuracy of about 0.1 mag can be ob- 
tained out to the Virgo and Fornax clusters of galaxies. A detailed comparison with well-known 
primary distance indicators is now required in order to establish the reliability and applicability 
of the WLR. In this context we have set to apply the methods of quantitative spectroscopy to A 
and B supergiants in galaxies with known, mostly Cepheid-based, distances. Among the galaxies 
studied within the HST Extragalactic Distance Scale Key Project (Kennicutt, Freedman & Mould 
1995, Mould et al. 2000), NGC 3621 has been chosen for its moderate distance (40% that of the 
Virgo cluster), maximizing our chances of detecting blue supergiants. In this Letter we discuss 
the identification of suitable targets and the testing of our spectroscopic tools with the current 
observational capabilities of an 8m telescope. 



For the spectroscopic target selection we obtained broad-band B, V and / CCD frames of 
NGC 3621 with Antu, the first unit of the Very Large Telescope (VLT) operated by the European 
Southern Observatory on Cerro Paranal, equipped with FORS (FOcal Reducer and low dispersion 
Spectrograph) in imaging mode. The crowded-field photometry was carried out with the IRAF 
version of daophot (Davis 1994), and indicates that the brightest stellar objects in this galaxy 



M Voo (R/R Q f 5 = L 



l/a 



(1) 



2. Target selection and observations 
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appear at V = 20, corresponding approximately to an absolute magnitude My = —10, given 
the Cepheid distance modulus (m — M) = 29.1 (±0.2) and average reddening E(V — I) = 0.3 
(Rawson et al. 1997). We selected isolated stars brighter than V = 22 and with color index —0.1 < 
(V — I) < 0.6, corresponding to the expected location of reddened blue supergiants in the colour- 
magnitude diagram. The target selection was finalized after inspection of archival HST WFPC2 
images of a smaller field in NGC 3621 (to avoid, where possible, contamination by neighbouring 
stars in our spectra), and rejection of objects showing nebular emission on top of stellar features, 
as revealed by short exposure spectra obtained during the commissioning phase of FORS2 at the 
second Unit Telescope of the VLT. Two H II regions (one of which located within a slit together with 
a stellar object) were included in the final sample, for future comparison of chemical abundances 
measured independently from stellar and nebular lines. 

Optical spectra of our candidates were obtained on March 1-2, 2000, under 0.8 arcsec seeing 
conditions. The FORS instrument possesses excellent multiplexing capabilities, allowing spectra of 

19 different faint objects to be taken simultaneously through separate slitlets, each approximately 

20 arcsec in length. Fig. 1 shows the distribution of the targets, and Table 1 summarizes their 
position and photometric data. Several exposures, adding up to a total integration time of 10.7 
hours, were obtained with a single multislit setup (600 gr/mm grism, 1 arcsec slits), covering as a 
minimum the 3,800-4,900 A wavelength range at 5 A resolution. This range contains most of the 
spectral features useful for deriving the main stellar parameters (effective temperature, gravity and 
chemical composition). 

The photometric selection proved to be rather successful. We have identified in our spectra 
ten supergiants, subdivided into spectral types as follows: B (one star), A (six), F (one), and 
two Luminous Blue Variable (LBV) candidates. Of the remaining targets, four have composite 
spectra, while four show nebular emission superposed on stellar features, despite our efforts to 
avoid such occurrences. In Fig. 2 one early A-supergiant (a), one early F-supergiant (b) and one 
LBV candidate (c) are displayed. The quality of the extracted spectra (with a S/N ratio of order 
50) is excellent, making the classification of spectral types using galactic standard stars, together 
with a first estimate of effective temperatures, straightforward. 

3. Stellar parameters and metallicity: a case study 

Using line-blanketed model atmospheres in hydrostatic equilibrium (Kurucz 1991; Kurucz 1992) 
and detailed NLTE/LTE line formation calculations (Przybilla et al. 2000a; Przybilla et al. 2000b) 
a quantitative analysis of the spectra is possible. We concentrate here on the object at the top in 
Fig. 2 (hereafter star a), for which we estimate an effective temperature T c gf = 9000±400 K from 
the spectral classification and subsequently derive a gravity log g = 1.05=L0.15 (cgs) from fitting the 
higher Balmer series members (from H7 upward). A microturbulent velocity typical for luminous 
supergiants, £ = 8 km s _1 , is chosen for the spectrum synthesis calculations. A more detailed anal- 
ysis comprising the remaining objects will be published elsewhere . The spectrum is well fit with a 
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chemical composition comparable to that of the Large Magellanic Cloud assuming a reduction of 
metallicity by a factor of two with respect to solar composition (Haser et al. 1998, see Fig. 3), while 
a significantly lower metallicity, as found in the Small Magellanic Cloud (0.2 solar), can confidently 
be ruled out. Previous work using the emission lines of H II regions close to our stellar target has 
found a reduction of the oxygen abundance of the same order (Ryder 1995). Our analysis provides 
additional information concerning elements such as magnesium, iron, chromium and titanium. 

It is important to note that the spectral resolution of 5 A allows a reasonable estimate of 
stellar metallicity by synthesizing the entire observed spectrum. We have tested this technique 
by reproducing the results obtained by detailed quantitative analyses of high resolution, high S/N 
spectra of Local Group A-supergiants (Przybilla, 2001, in prep.). However, we stress that it is 
difficult to determine precise effective temperatures from ionization equilibria of weak lines and 
microturbulence velocities and abundances from individual lines. Here we have to rely on the 
experience obtained from the analyses of Local Group objects using higher resolution spectra. In 
the next phase we intend to work at a resolution of 2 A as provided, for instance, by FORS 2 at 
the VLT. 

We determined the interstellar reddening by fitting the calculated energy distribution to the 
fluxes corresponding to the measured magnitudes. A total E(B— V) = 0.12 is found, consistent with 
most, if not all, of the effect as being due to galactic foreground reddening (Burstein & Heiles 1982). 
This value is smaller than the one obtained by the HST Key Project {E(B — V) = 0.23, Rawson et 
al. 1997), however their WFPC2 field is closer to the central part of the galaxy, where dust lanes 
are prominent (Fig. 1). We therefore demonstrate that, in principle, stellar spectroscopy allows 
us to map interstellar extinction for accurate distance determinations. From the Cepheid distance 
modulus and an interstellar extinction Ay = 0.37 we derive an absolute magnitude My = —9.0 
(±0.2) for our object. The corresponding stellar radius, calculated from the model atmosphere flux 
in the V filter, is R = 250 ±25 R Q , yielding a luminosity L = 3.7 x 10 5 L Q . 

4. Wind parameters and the WLR 

Star a, with its very high luminosity, is the ideal target for testing whether it is possible to 
measure wind momenta of A-supergiants in galaxies far beyond the Local Group. Mass-loss rates 
M and terminal velocities Voo are usually determined from the shape of the Ha stellar wind line 
profile. However, because of the enormous strength of the wind driven by the luminosity, also 
H/3 is affected by mass-loss, which allows a first estimate of the stellar wind momentum. As in 
previous work in the Milky Way (Kudritzki et al. 1999) and M31 (McCarthy et al. 1997) we adopt 
"unified model atmospheres", which allow for departures from Local Thermodynamic Equilibrium, 
are spherically extended and include the hydrodynamic effects of stellar winds in their stratification 
(Santolaya-Rey, Puis & Herrero 1997) . The calculated H/3 profiles in Fig. 4 show very clearly how 
mass-loss affects the profile shape. The absorption profile is significantly filled by wind emission, 
enabling a reasonable determination of M. Since, contrary to Ha, the H/3 fit does not allow the 
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direct determination of terminal velocities, we adopt a typical value Voo = 200 km/s. Comparing 
with other A-supergiants (Kudritzki et al. 1999; McCarthy et al. 1997) the uncertainty introduced 
by this assumption is of the order of 20%. The mass- loss of the A-supergiant obtained in this way 
is 3 x 10 -6 M & yr _1 , while the corresponding modified stellar wind momentum flow is 5.9 x 10 28 
erg cm -1 , with an accuracy of roughly 30%. 

Fig. 4 compares with the WLR for galactic and M31 supergiants of solar metallicity, for which 
we have adopted the stellar parameters given by Venn et al. (2000), Kudritzki et al. (1999) and 
McCarthy et al. (1997), together with a more recent distance determination to M31 equal to 
783 Kpc (Holland 1998), and a revised photometry for the M31 supergiants (Venn et al. 2000a). 
Because of the small bolometric corrections of A-supergiants we use absolute magnitude rather 
than luminosity to display the relationship. Star a is located slightly below the regression line. 
A natural explanation for this is the somewhat lower metallicity inferred from the analysis of 
the photospheric spectrum, which leads to lower wind momenta (Kudritzki 1998; Kudritzki 2000; 
Puis et al. 1996; McCarthy et al. 1995; Haser et al. 1998). An empirical calibration of the WLR 
metallicity dependence using A-supergiants in the Magellanic Clouds is presently under way. The 
slight discrepancy might also be due to the uncertainty of the H/3 fit. The forthcoming VLT 
observations of all of our targets at Ha wavelengths will allow us to answer this question and to 
put additional data points onto the WLR-plane. We are confident that this will also provide an 
independent constraint on the distance of NGC 3621. 

We thank the many people within the FORS consortium, especially B. Muschielok and I. 
Appenzeller, for the FORS2 commissioning data used in the preliminary phases of this work. 
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Table 1. FORS spectroscopic targets 



Slit no. 


R.A 


. (J2000) 


DEC (J2000) 


V 


B — V 


V — I 


Comments 


i 
i 


1 1 

J. 1 


18 


14.3 


-32 


45 


40.1 




n i a 


n n 

u.oo 




2 


11 


18 


20.6 


-32 


45 


31.4 


20.41 


0.12 


0.28 


LBV candidate 


3 


11 


18 


10.3 


-32 


46 


25.2 


21.67 


0.00 


0.17 




4 


11 


18 


27.2 


-32 


46 


02.6 


20.43 


0.32 


0.43 


star b 


5 


11 


18 


15.3 


-32 


46 


54.1 


20.76 


0.07 


0.34 




6 


11 


18 


07.9 


-32 


47 


45.9 


21.80 


0.11 


0.16 




7 


11 


18 


15.8 


-32 


47 


39.2 


21.06 


0.28 


0.28 




8 


11 


18 


15.9 


-32 


48 


11.2 


20.19 


0.30 


0.40 




9 


11 


18 


28.7 


-32 


48 


04.6 


20.47 


0.17 


0.28 


star a 


10 


11 


18 


10.7 


-32 


49 


11.9 


22.32 


0.03 


0.05 




11 


11 


18 


30.5 


-32 


48 


37.9 


21.87 


-0.02 


0.01 


H II region 


12 


11 


18, 


18.1 


-32 


49 


31.8 


20.39 


0.25 


0.48 


star c, LBV candidate 


13 


11 


18 


16.7 


-32 


50 


13.1 


21.07 


0.45 


0.60 


H II region + star 


14 


11 


18 


26.6 


-32 


50 


05.8 


21.90 


0.20 


0.18 




15 


11 


18 


20.2 


-32 


50 


35.2 


21.68 


0.15 


0.31 




16 


11 


18 


23.3 


-32 


51 


02.7 


21.74 


0.07 


0.05 




17 


11 


18 


27.5 


-32 


51 


05.4 


21.01 


0.21 


0.33 




18 


11 


18 


18.2 


-32 


51 


58.9 


21.04 


0.00 


0.09 




19 


11 


18 


19.6 


-32 


52 


10.1 


19.51 


0.79 


1.05 
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Fig. 1— The 19 objects observed with the VLT and the FORS spectrograph in NGC 3621 arc 
marked with circles (stars) and squares (H II regions) on a color image obtained by combining 
5-minute B-, V— and /—band frames taken with the same instrument in imaging mode. The field 
of view is approximately 7x7 arcmin. Different coloured markers are used for clarity only. The 
stars labeled as a, b and c are further described in Fig 2. 



Fig. 2. — Examples of continuum-normalized stellar spectra obtained in NGC 3621, displayed along 
with galactic supergiant stars of similar spectral morphology and at the same resolution. Ordinate 
ticks are spaced by 0.5 continuum units, (a) An early A-supergiant (red, star a in the text), 
bracketed by galactic AOIa and A2Ia template spectra (black). Line identification for the strongest 
features is shown at the top (unlabeled shorter marks are used for Fe II lines), (b) An early F- 
supergiant (red) compared to galactic FOIa and F2Ia template spectra (black), (c) A Luminous 
Blue Variable (LBV) candidate. Notice the prominence of several Fe II emission features (identified 
by the vertical marks), and the presence of nebular emission, indicated by the hydrogen Balmer 
lines in emission. 



Fig. 3. — The observed, continuum-normalized spectrum of the A-supergiant star a (in black) in the 
4270-4650 A wavelength range is compared to the model predictions (T e g = 9000 K, log g = 1.05) 
at LMC (red) and SMC (blue dotted) chemical abundances, corresponding to 0.5 and 0.2 solar, 
respectively. Line identification is provided (Fe II lines are indicated by the shorter marks). 



Fig. 4. — The relationship between wind momentum (in cgs units) and luminosity (expressed in 
terms of absolute V magnitude) for A-supergiants in the Galaxy (open circles), M31 (filled circles) 
and star a (shown here by the error bars). The latter was excluded in the calculation of the 
regression line. The inset shows how the mass-loss rate was estimated by means of H/3 line profile 
fits with NLTE unified wind model atmospheres including stellar winds. The vertical and horizontal 
bars represent 0.05 continuum units and 5 A, respectively. The observed spectrum is shown by the 
blue line. Our final profile fit, corresponding to a mass-loss rate of 3 x 10 _6 M Q yr" 1 , is shown in 
red. Increasing the mass-loss rate to 4 x 1O _6 M yr _1 leads to the upper dotted line, showing an 
incipient P-Cygni profile. Lower mass-loss rates (2 x 10~ 6 and 1O~ 7 M yr _1 ) lead to the additional 
dotted profiles. At least part of the discrepancies between the observed and theoretical profiles can 
be attributed to the presence of metal lines (at wavelengths indicated by the vertical marks) not 
accounted for in the stellar wind line profile calculations. The deviations are however within the 
observational uncertainties. 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0012372vl 
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